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ABS T RACT. 

The CAP-TSD (Computational Acroelasticily 
Program - Transonic Small Disturbance) code, developed 
at the NASA - Langley Research Center, is applied to the 
Active Flexible Wing (AFW) wind-tunnel model for 
prediction of the model's transonic aeroclastic behavior. 
Static aeroelastic solutions using CAP-TSD are 
computed. Dynamic (flutter) analyses are then performed 
as perturbations about the static aeroclastic deformations 
of the AFW. The accuracy of the static aeroclastic 
procedure is investigated by comparing analytical results 
to those from previous AFW wind-tunnel experiments. 
Dynamic results are presented in the form of root loci at 
different Mach numbers for a heavy gas and air. The 
resultant flutter boundaries for both gases arc also 
presented. The effects of viscous damping and angle-of- 
attack, on the flutter boundary in air, are presented as 
well. 

INTRODUCTION 

An understanding of the aeroclastic behavior of flight 
vehicles in the transonic regime is of great importance for 
flight safety. For example, it is well known that aircraft 
flying into or through the transonic regime may encounter 
a region of reduced flutter speed known as the transonic 
flutter dip. Valuable insight into the nature of this 
transonic flutter dip phenomena is provided by Isogai 1 for 
a two-dimensional airfoil, while comparison of 
aerodynamic theory with the experiments reported by 
Davis and Malcolm^ reveal the limitations of linear 
theory applied in the transonic regime. Linear 
aerodynamics, although highly successful in the subsonic 
and supersonic regimes, cannot normally be used to 
accurately predict transonic aeroclastic behavior. 
Transonic flow equations capable of modelling flow 
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nonlincaritics (shocks, boundary layer, separation and 
vorticiiy) and boundary condition nonlinearities (airfoil 
thickness and shape, and large deflections) must then 
be solved. The surveys by Edwards and Thomas^ and 
Ballhaus and Bridgcman 4 review recent computational 
developments in the field of transonic aeroelasticity. 
Some of these developments include modelling of the 
Navier-Stokes equations 5 and the Euler equations 6 for 
flutter analysis. Application of these higher order 
formulations, however, has primarily been limited to two- 
dimensional configurations, due to the large 
computational costs incurred. Certain assumptions 
regarding the flow can be made to yield reduced order 
formulations such as the full-potential equation 7 and the 
computationally efficient transonic small-disturbance 
(TSD) equation. Research efforts involving the TSD 
formulation include the development of the XTRAN3S 
code**, the work by Yang, Guruswamy, and Striz^, and 
many others. 

A transonic aerodynamics code known as CAP-TSD 
(Computational Aeroelasticity Erogram-Transonic Small 
Disturbance) has been developed at the NASA - Langley 
Research Center (LaRC). CAP-TSD is capable of 
handling multiple lifting surfaces with control surfaces, 
bodies (nacelles, pylons, stores), vertical surfaces, and a 
fuselage, and solves the TSD equation using an efficient 
approximate factorization scheme*®. References 11-12 
verified the code's ability to accurately predict steady and 
unsteady pressures for wings and configurations at 
subsonic, transonic, and supersonic Mach numbers. 
Flutter prediction using CAP-TSD for two thin, swept- 
and-tapered wings compared well with experimental flutter 
results* 3. The goal of the present study was to define the 
transonic flutter boundary of the Active Flexible Wing 
(AFW) wind-tunnel model 14 * 15 , for use as guidance 
during flutter testing, and to evaluate CAP-TSD’s flutter 
prediction capability for a complete and realistic aircraft 
configuration. 

The Active Flexible Wing (Fig. 1) model is a full- 
span, sting-mounted wind-tunnel model designed and built 
by the Rockwell International Corporation. The main 
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Fig. 1 The AFW in NASA-LaRC’s Transonic 
Dynamics Tunnel (TDT). 


goal of the present AFW project is to design, implement 
and validate digital control laws for flutter suppression ^ 
and roll maneuver load alleviation. A priori knowledge of 
possible regions of instability are, therefore, crucial. 

This paper first presents the computational procedures 
incorporated in CAP-TSD. This includes a brief 
description of the TSD formulation and the coupled 
aerodynamic and structural equations of motion that are 
integrated in time. These equations are used for both 
static aeroelastic and dynamic analyses of the AFW. An 
important conclusion of the studies by Yates, Wynne, and 
Farmer 12 and Yates and Chu^ was that the accuracy of 
the transonic flutter prediction is highly dependent on the 
accuracy of the static aeroelastic state of the wing. As a 
result, a procedure for computing static aeroelastic 
deformations is presented. The dynamic behavior is 
computed as a perturbation about previously computed 
static aeroelastic solutions. The resultant dynamic time 
histories of the generalized displacements are then 
analyzed using a modal identification technique to 
estimate the stability parameters (root locus) of the 
system at a given Mach number and dynamic pressure. 
Dynamic results arc presented in the form of root locus 
excursions at different Mach numbers in a heavy gas 
(Freon* 12) and in air. Flutter boundaries for the heavy 
gas and air, variations in angle-of-altack and viscous 
damping, and comparisons with experimental flutter 
results are also presented. 

COMPUTATIONAL PROCEDURES 
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small-disturbance (TSD) equation. The TSD potential 
equation is defined by 

Mi «> t + 24» x ) t = [(1 - Mi + F<D l + Gc>y ] x + 

(0y + H^y )y + ) Z (1) 


where M^ is the Mach number and $ is the disturbance 
velocity potential. 

Several choices are available for the coefficients F, 
G, and H depending upon the assumptions used in 
deriving the TSD equation. For transonic applications, 
the coefficients are herein defined as 

F=-y(Y+l)Mi . 

G = l(Y-3)Mi, 

H = - (Y - 1) Mi (2) 

The linear potential equation is recovered by simply 
setting F, G, and H equal to zero. 

Equation (1) is solved within CAP-TSD by a time- 
accurate approximate factorization (AF) algorithm 
developed by Batina 10 . In Refs. 10 to 12, the AF 
algorithm was shown to be e ffic ient for appli cation lo 
steady or unsteady transonic flow proWems,.. It can 
provide accurate solutions in only several hundred time 
steps yielding a significant computational cost savings 
when compared to alternative methods. Several algorithm 
modifications have been made which improve the 
stability of the AF algorithm and "the accuracy of the 
results 19 * 20 . The CAP-TSD program can treat 
configurations with combinations of lifting surfaces and 
bodies including canard, wing, taif control surfaces, tip 
launchers, pylons, fuselage, stores, and nacelles. 

The configuration capability of ihe current version of 
CAP-TSD permits the calculation of pressures on the 
fuselage and bodies. In the present study, modal 
perturbations of the fuselage and bodies arc not included 
in the boundary conditions and the integration of the 
pressures for the generalized aerodynamic forces of the 
fuselage and bodies are not included in the a eroelas tic 
solution. However, the aerodynamic influence of both 
the fuselage and wing tip body of the AFW model arc 
included as interference effects upon the wing pressures. 


In this section, an overview of the computational 
procedures is presented including a description of the 
CAP-TSD program, the aeroelastic equations of motion, 
the time-marching solution of these equations, and the 
modal identification of the resulting free decay transients. 

CAP-TSD Program 


Equations of Motion 

The aeroelastic equations of motion are based on a 
right-hand orthogonal coordinate system with the x- 
direction defined as positive downstream, y-direction 
positive out the right wing, and the z-dircction positive 
upward. The equations of motion may be written as 


Mq + Cq + Kq = Q 


(3) 


The CAP-TSD program is a finite-difference program 
which solves the general-frequency modified transonic 
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where q is a vector of generalized displacements, M is the 
generalized mass matrix, C is the damping matrix, and K 
is the stiffness matrix. Q is the vector of generalized 
forces where its elements are defined by 


Qi = 


pu 2 , f A P h i dS 

2 J 2 2 

1 \ P tr/2 c r 


and Ap is the lifting pressure, p is the fluid density, c r is 
the root chord, U is die freestream velocity, S is the area 
of the lifting surface(s) and hj is the vibration mode 

shape. Equation (3) is rewritten as 

\ 

q = -M' 1 Kq - M' 1 Cq + M' 1 Q (4) 

to permit integration of the equation with respect to time. 

Time-Marching Aeroelastic Solution 


The aeroelastic solution procedure implemented 
within CAP-TSD for integrating Eq. (4) is similar to that 
described by Edwards, Bennett, Whitlow, and Seidel^ 
Equation (4) is composed of normal mode equations 
which may be expressed in linear, first-order state-space 
form as 


x. = Ax. + Bu 

i i i 


where 

x i = K \ i T 


and 


A = 


0 1 ~ 
-m' 1 k. -m' 1 c. 

ii ii 


(5) 


B = m' 1 

i 



0 

1 


J AC p h. dS/c? 


AC p 


Ap 

p U 2 /2 


In these definitions, mj, Cj, and kj arc elements of the 
mass, damping, and stiffness matrices, respectively, 
corresponding to mode i. The analytical solution to Eq. 
(5) and a description of its numerical implementation in 
CAP-TSD is found in Refs. 13 and 21. 

For aeroelastic analysis, two steps arc generally 
required in performing the calculations. In the first step, 
the steady-state flow field is calculated to account for wing 


thickness, camber, mean angle-of-attack, and static 
aeroelastic deformation, thus providing the starting flow 
field for the dynamic aeroelastic analysis. Previously 
published CAP-TSD flutter studies analyzed only 
symmetric airfoils at zero angle-of-attack^, thereby 
avoiding the problem of static aeroelastic deformations. 
For the AFW, the wing is unsymmetric and rigged at a 
non-zero angle-of-attack, so a procedure for computing 
static aeroelastic solutions had to be developed before an 
accurate dynamic analysis could be performed. The 
dynamic analysis would then be a perturbation about a 
converged static aeroelastic solution at each Mach number 
and dynamic pressure of interest. 

The procedure developed and applied in this study for 
computing static aeroelastic deformations is to allow the 
structure and aerodynamics to interact with no initial 
excitation (no initial deflection or velocity) and with a 
large value of viscous damping to prevent divergence of 
the solution. This method resulted in convergence of the 
generalized displacements. Static aeroelastic deformations 
should be independent of viscous damping and different 
values of viscous damping (£=.375, .707, and .99) were 
evaluated. A typical result for this type of analysis is 
presented in Fig. 2, which shows a representative 
variation of a generalized displacement as a function of 
computational time steps for three values of viscous 
damping. It is clear from Fig. 2 that the convergence is 
indeed independent of the value of viscous damping. 



Fig. 2 Convergence of generalized displacements 
for different values of viscous damping. 


Furthermore, the larger the value of viscous damping, 
the faster the convergence. Therefore, the highest value 
of viscous damping (£ = 0.99) was used in order to 
accelerate the static aeroelastic solution. For the 
applications presented herein, 1000- 2000 time steps were 
used to converge the static aeroelastic solutions. An 
interesting result of this procedure is that it allows the 
computation of static aeroelastic deformations at dynamic 
pressures above the fluuer dynamic pressure for the AFW. 

Once converged static aeroelastic solutions are 
computed, the next step is to prescribe an initial 
disturbance to begin the dynamic structural integration. 
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Disturbance (or modal) velocities in the first three modes 
are used as initial perturbations. About 7 cycles of the 
lowest frequency (first) mode were needed for accurate 
modal identification. For a constant, non-dimensional 
time step of .01 , this required 4000 time steps in the 
heavy gas and 8000 time steps in air. In determining a 
flutter point, the freestream Mach number. Moo, and the 
associated freestream speed, U, were held fixed. A value 
of the dynamic pressure pU^/2 is then used and free decay 
transients are computed. These resulting transients of the 
generalized coordinates are analyzed for their content of 
damped or growing sine-waves, with the rates of growth 
or decay indicating whether the dynamic pressure is above 
or below the flutter value. This analysis then indicates 
whether to increase or decrease the value of dynamic 
pressure in subsequent runs to determine a neutrally stable 
result. 




b) identified modal components 
and offset 

Fig. 3 Example of dynamic decay record and its 
modal components. 

Modal Id entification 

As previously mentioned, CAP-TSD generates 
free decay transients that must be analyzed for the modal 
stability characteristics. A typical transient for the AFW 
model, calculated using CAP-TSD is shown in Fig. 3(a). 
The first three modes used in the analysis were excited by 
specifying an initial condition for each modal velocity to 
produce a complex decay record. This record is analyzed 
using a least-squares curve-fit of the response data with 
complex exponential functions. The program utilized is a 
derivative of the one described in Ref. 22. The 
components of the transient of Fig. 3(a) are plotted in 
Fig. 3(b) to the same scale. The free decay properties of 
each mode for this condition are readily apparent and the 
mean or offset value is the static acroelaslic deformation 
of the mode being analyzed. A sufficient range of 


dynamic pressure must be considered to determine all 
relevant flutter points. 

AMALYSIS AND RESULTS 

CAP-TSD ComputatitmaJL. Model 

The AFW geometry data was obtained from Rockwell 
International, including detailed airfoil shape information. 
From this geometry data, a half-span model, with 
symmetry specified at the centerline, was generated. This 
CAP-TSD model consists of a fuselage, the addition of 
the region aft of the main wing and next to the fuselage 
referred to as the coat-tail, the main wing with all four 
control surfaces, and the wing tip ballast store. The grid 
dimensions for this model are 134x51x62 in the x-, y-, 
and z-directions respectively for a total of 423,708 grid 
points. The grid extends 10 root chords upstream, 10 root 
chords downstream, 2 semi-span lengths in the y- 
direction, and 10 root chords in the positive and negative 
z-direction. Modelling of the wind-tunnel sting mount is 
done by extending the computational fuselage aft to the 
downstream boundary. The grid density is increased in 
regions where large changes in the flow are expected, such 
as at the leading edge, trailing edge, wing tip, and 
control-surface sides and hinge lines. The four control 
surfaces are the leading-edge inboard (LEI), leading-effje 
outboard (LEO), trailing-edge inboard (TEI), and trailing- 
edge outboard (TEO). Each control surface has a chord 
that is 25% of the local chord and a span that is 28% of 
the semi-span. The airfoil definition includes the control 
surface actuator bumps on the outboard half of the wing. 
There also exist slight discontinuities on the wind-tunnel 
model where wing box and control surfaces meet (at the 
quarter- and three-quarter chord). These discontinuities are 
not included in the analytical model because of potential 
numerical difficulties. The effect of the actuator bumps 
and the control surface/wing box discontinuities on the 
measured and computed static pressure distributions will 
be presented in a subsequent section. A computer- 
generated picture of the CAP-TSD model of the AFW is 
shown in Fig. 4. Although not shown in the figure, a 



Fig. 4 CAP-TSD computational model of the AFW. 
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protrusion on the underside of the fuselage that houses the 
model's pitch actuator is also included in the analytical 
model. 

Analytical modes and frequencies were obtained from a 
finite-element model and separated into symmetric and 
antisymmetric modal data sets. The flutter analysis was 
performed using analytical mode shapes with measured 
frequencies (ground vibration test). The symmetric data 
was shown by linear analysis^ to be the most flutter 
critical in the higher, subsonic Mach number regime and 
so only these were analyzed in the present study. A total 
of eight symmetric modes were included in the model. 
The interpolation of mode shape displacements and slopes 
at the computational grid points is done via a surface 
spline^. Each structural section was splined separately 
and then recombined to form the necessary input to CAP- 
TSD. The separate structural sections are the wing box, 
coat-tail, and the four control surfaces. Slender bodies 
such as the fuselage and tip ballast store are not given any 
modal definition in CAP-TSD, as was previously 
mentioned, therefore no modal data was needed for these 
components. 


Sialic Aeradastic J ksiil t s 

The accuracy of the static aeroelastic solution can 
be determined by comparing analytical results with 
existing experimental data. There were three sets of 
experimental data, from the previous AFW tests in the 
heavy gas, available for this purpose. This data included : 
1) pressure coefficient distributions; 2) control-surface 
effectiveness parameters; and 3) static deflection data 
computed from experimental pressure distributions due to 
control-surface deflections^ 4 . It should be noted that the 
comparisons with the second and third sets of 
experimental data are not a direct assessment of the static 
aeroelastic procedure alone, since the accuracy of the 
control surface modelling within CAP-TSD is obviously 
an integral part of the result. Static deformation data with 
no control surface deflection is desirable, but, 
unfortunately, not readily available. Both sets of data, 
however, arc useful in observing the trends and behavior 
of the static aeroelastic procedure as well as the control 
surface modelling within CAP-TSD. Note that the AFW 
configuration for these previous tests did not include the 
lip ballast store used in the recent test. For the CAP- 
TSD calculations to compare with the earlier experiments, 
the tip ballast store was deleted and the tip fairing added. 

Pressure distributions- Figure 5 presents pressure 
coefficient distributions versus percent chord for CAP- 
TSD and experiment at Moo = 0.9 and a dynamic pressure, 

q, of 150 psf at the three span stations shown, where q is 
the percent semi-span. 

The overall agreement between analysis and 
experiment is good, with some discrepancies occuring 
near the trailing edge and wing tip. The first two span 
stations compare remarkably well from the leading edge 
up to about sixty percent of the local chord. Sudden 
changes in the flow can be seen near the quarter-chord at 


the second span station and near the three-quarter chord for 
all three span stations. These disruptions in the flow may 
be caused by the previously-mentioned physical 
discontinuities where wing box and control surfaces meet. 
At the second and third span stations, the effect of the 
actuator bumps on the lower surface pressures is evident. 
Agreement between analysis and experiment deteriorates at 
the third span station, possibly due to separated and/or tip 
vortex flow around the wing tip region. 

Analytical and experimental pressure data were 
also compared at a lower dynamic pressure (q=36 psf) 
although not presented herein. Since the static aeroelastic 
deformations at the higher dynamic pressure (q=150 psf) 
are larger than at the lower dynamic pressure (q=36 psf), 
the results at the higher dynamic pressure (Fig. 5) provide 
a more stringent test of the sialic aeroelastic procedure. 
Surprisingly, there exists an overall improvement in 
correlation between analytical and experimental pressure 
distributions at the higher dynamic pressure. This is 
possibly due to the fact that at the higher dynamic 
pressure, the flow tends to remain attached over a larger 
portion of the wing than at the lower dynamic pressure, 
creating a condition that is closer to the TSD assumptions 
of inviscid, attached flow. 



Fig. 5 Comparison of pressure distributions at Mc*=0.9 
and q=150 psf in the heavy gas. 
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At Moo = 0.95, q=36 psf, the analytical and measured 
pressure distributions differ significantly (not shown 
here), specifically, in the shock strength and location. 
Typical for isentropic, inviscid flow theory, the shock is 
predicted too far aft and too strong when compared with 
experimental results. Even the use of vorticity and 
entropy corrections in the analysis did not improve the 
results significantly. There arc evidently other nonlinear 
flow effects (separation, boundary layer) and a greater 
sensitivity to differences in the computational and 
physical airfoil shape that weaken the comparison at this 
test condition. 

Control-surface effectiveness - For the control- 
surface effectiveness parameters, the present study 
investigates only the lift coefficient due to control-surface 
deflection for all four control surfaces at one Mach 
number and dynamic pressure. T his was not meant to be 
an extensive investigation into the control surface 
capabilities^ of CAP-TSD as a detailed analysis of this 
capability would require a grid convergence study to 
determine the effects of varying grid densities in the 
vicinity of the control surfaces’ sides and hinge lines. 
Furthermore, with deflected controls, boundary layer 
effects are expected to be significant and affect the 
effectiveness of the controls whereas the current 
calculations are far inviscid, attached flow. 

The analysis was performed in much the same way 
that the data was measured in the wind tunnel. For each 
control surface deflection, the model was allowed to 
converge to a static aeroclastic solution and the resulting 


lift coefficient was recorded. Experimental data was 
available at +5, -5, and 0 degrees of control surface 
deflection. Figure 6 presents the comparison between 
CAP-TSD and experimental control surface effectiveness 
parameters for all four control surfaces at Moo = 0.9 and 
q= 150 psf. 

It can be seen that the trailing-edge controls arc 
analytically more effective than the experimental results, 
while the leading-edge controls are analytically less 
effective. It is quite possible that at the leading-edge, 
vortex flows induced by the control surface deflections are 
enhancing the experimental effectiveness of the control 
surfaces. At the trailing edge, however, separated and 
boundary layer flows may be reducing the measured 
control surface effectiveness, most noticeably outboard. 
The fact that the LEI control surface correlates the best 
with experiment at this flight condition is consistent with 
the fact that the best pressure correlations are at the 
leading edge and inboard of the wing. The TEI control 
surface results are surprisingly good, mostly due to the 
large loads generated by this surface, thereby reducing the 
relative error. Comparison between analysis and 
experiment for the TEO control surface is poor, probably 
due to viscous (boundary layer) and tip vortex flows, not 
modelled in TSD theory. 

In general, the prediction of control surface 
effectiveness parameters by CAP-TSD can qualitatively 
identify the most effective (TEI) and least effective (TEO) 
control surface, but improvement is needed for 
quantitative prediction. 



Inboard Outboard 



Fig. 6 Control surface effectiveness, Moo=0.9, 
q=150 psf in the heavy gas. 


Static Aeroelastic Deflections - The third and final 
set of data is static aeroelastic deflections. In order to 
determine wing deflections, a distributed, coarse grid 
of vertical forces using experimental pressure 
distributions was computed^. These forces were then 
multiplied by the model's measured structural flexibility 
matrix to obtain wing deflections. These deflections are 
therefore quasi-expcrimental, since they were not measured 
directly. The CAP-TSD wing deflections were computed 
by the summation of the products of each converged 
generalized displacement (from the static aeroelastic 
solution) with its corresponding mode shape. Figure 7 
gives the comparison of the quasi-experimental and the 
CAP-TSD wing deflections along the chord, or wing 
station, at three spanwise locations due to a TEO 
deflection of plus five degrees (trailing-edge down) at Moo 
= 0.9, q=l 50 psf. It should be mentioned that those arc 
the resultant modal deflections and not the actual shape of 
the wing, since the latter should include the downward 
deflection of the TEO control surface. 

The larger analytical deflections are consistent with 
the larger analytical loads predicted for the TEO control 
surface (Fig. 6). The nature of the wing deflection is well 
described with both analysis and experiment showing a 
bend-up, twist-down deformation. The corresponding 
results for a TEO control surface deflection of minus five 
degrees, not presented here, compared better in terms of 
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deflection, inches 



Fig. 7 Comparison of quasi-experimental and 
analytical static aeroeiastic deflections, 

Moo=0.9, q=150psf in the heavy gas. 

the magnitude of deflections and moderately well in 
predicting the amount of twist at the wing tip. However, 
the resulting deformations (quasi-experimental and 
calculated) were much smaller for the minus five degree 
control surface deflection. 

Based upon these results, the static aeroeiastic 
solution is viewed as a reasonably accurate approach. The 
accuracy of the solution, when combined with control 
surface deflections, is diminished although the trends 
remain within reason. Methods for improving the static 
aeroeiastic solution include the application of measured 
mode shapes (the present aeroeiastic analyses uses 
analytical mode shapes with measured frequencies) and an 
increased number of mode shapes. Improvement of the 
control surface modelling may be possible by a finer grid 
at the control surface boundaries and limiting application 
to small deflection angles. This may only be valid for the 
leading-edge controls, since incorporation of boundary 
layer and separated flows will probably be needed for 
improvement in the modelling of the trailing-edge control 
surfaces . 

Dynamic Results 

Results in the Heavy Gas- The root locus of the 
first four clastic modes with nonlinear aerodynamics for 
Moo = 0.9 and 1.5 degrees anglc-of-attack with no viscous 
damping is presented in Figure 8. It should be mentioned 
that the wing tip ballast store is included in these results. 
The flutter mechanism involves the coalescence of the 
second (first bending) and third (first torsion) elastic 
modes rcsluting in a flutter dynamic pressure of 213 psf 
and a flutter frequency of 9.7 Hz. The branch of the root 
locus for first bending yields the flutter condition. 

The resultant flutter boundary, Mach number versus 
flutter dynamic pressure, and corresponding flutter 
frequencies, are shown in figure 9. A moderate transonic 


frequency, rad/sec 



damping, 1/scc 


Fig. 8 Root locus of first four elastic modes with 
nonlinear aerodynamics, M©o=0.9, alpha= 

1 .5 deg, and no viscous damping. 

flutter ’’dip" is evident, with the bottom of the "dip" ai 
Moo = 0.95, a dynamic pressure of 153 psf and a flutter 
frequency of 9.2 Hz. 

It is interesting to note that, in the flutter analysis 
with doublet lattice (linear) aerodynamics^ , the predicted 
flutter mechanism is a torsion-dominated (third mode) 
instability for all Mach numbers analyzed in the heavy 
gas. The nonlinear aerodynamic terms are, therefore, 
altering the nature of the flutter mechanism from a 
torsion-dominated instability to a bending-dominated 
instability, as seen in Fig. 8. 



Fig. 9 Flutter boundary wiLh nonlinear aerodynamics 
at alpha=1.5 deg and no viscous damping in 
the heavy gas. 
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Fig. 10 Root locus of first four clastic modes with 
nonlinear aerodynamics at Moo=0.5, alpha= 
L5 deg and no viscous damping, in air. 

Results in Air - Figure 10 is the root locus competed 
using nonlinear aerodynamics at Moo = 0.5 and 1.5 degrees 
angle-of-attack with no viscous damping. Although the 
aerodynamics at = 0.5 are linear, the analysis was 
performed using the nonlinear aerodynamic equations so 
that the effect of the nonlinear terms on the flutter 
mechanism could be evaluated as Mach number was 
varied. Figure 10 shows the coalescence between the 
second and third modes, with the third mode (torsion) 
dominating the instability at a dynamic pressure of 245 
psf and a flutter frequency of 11.14 Hz. The flutter 
analysis using linear aerodynamics in air *5 also predicted 
a torsion-dominated instability for all Mach numbers from 
Moo = 0.5 to Moo = 0.9 (highest Mach number evaluated). 
This implies that the nonlinear terms have little effect on 
the aerodynamics at Moo = 0.5, as expected. 

The root locus for Moo = 0.85 with nonlinear 
aerodynamics at 1.5 degrees angle-of-attack and no viscous 
damping is shown in Fig. 11, with a flutter dynamic 
pressure of 204 psf and a flutter frequency of 9.55 Hz. At 
Moo = 0.85, a reversal of the dominant nutter mode from a 
torsion-dominated instability (Moo = 0.5) to a bending- 
dominated instability is noticed. Since the dominant 
flutter mode, predicted with linear aerodynamics, does not 
vary with Mach number, this reversal of the dominant 
flutter mode between Moo = 0,5 and Moo = 0.85 (and at 
M OO — 0.9 in the heavy gas), is a measure of the 
increasing sensitivity to differences in the aerodynamic 
modelling. As Mach number is further increased, the 
dominance of the bending flutter mode continues to grow, 
as if approaching a singlc-dcgrcc-of-frccdom (bending) 
instability. This is consistent with the result by Isogai^. 
Reference 1 shows that for a two-dimensional wing (with 
vibrational properties similar to those of a typical, 
streamwise section of an aft-swept wing), the flutter 
mechanism at subsonic Mach numbers is the classical 


bending-torsion instability, driven by the phase lag 
between the two modes. As Mach number is increased, 
however, the phase lag between the two modes is 
gradually reduced to near zero, signaling the presence of a 
single-degrec-of-freedom motion, or bending in the case of 
Ref. 1. 

Figure 12 is a plot of Mach number versus flutter 
dynamic pressure at 1.5 degrees angle-of-attack and no 
viscous damping. A severe transonic flutter "dip" is 
present, with the bottom of the "dip" at Moo = 0.93, a 
dynamic pressure of 20 psf and a flutter frequency of 7.16 
Hz. A variation in transonic flutter "dip" between air and 
the heavy gas is expected since, as Mykytow^ pointed 
out, the greater the mass ratio, the greater the transonic 
flutter "dip". Reference 1 supports this statement by 
showing that increases in mass ratio drive the aeroelastic 
system towards an earlier onset of the nearly singlc- 
degree-of-freedom instabiliy. The AFW model in air 
experiences about double the mass ratio experienced in the 
heavy gas, and the effect of this increase in mass ratio can 
be seen in Figs. 8 and 11. The Mo© = 0.85 flutter 
mechanism in air (Fig. 11) exhibits a slightly stronger 
bending-dominated instability than the Moo = 0.90 flutter 
mechanism in the heavy gas (Fig. 8). However, the 
magnitude and steepness of the "dip" in air is surprising. 

Damping and ang le-of-attack variations- The 
structural damping of the AFW model was determined 
from GVT tests to be about 1.5 % critical damping. In 
order to account for this, flutter analyses were performed 
at Mo© = 0.5, 0.9, 0.93, and 0.95 with a viscous damping 
value of 0.015 for each mode. Note that the static 
aeroelastic analysis for these cases did not have to be 

frequency, rad/sec 



damping, 1/sec 

Fig. 1 1 Root locus for the first four clastic modes 
with nonlinear aerodynamics, Moo=0.85, 
alpha=1.5 deg and no viscous damping in air. 
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Fig. 12 Effect of viscous damping and angle-of- 
attack on the flutter boundary in air. 


rerun, since the viscous damping affects only the dynamic 
analysis. The flutter boundary including a viscous 
damping of 0.015 at 1.5 degrees angle-of-attack can also 
be seen in Fig. 12. The bottom of the "dip" is higher but 
still at Moo = 0.93, going from a dynamic pressure of 20 
psf with no damping to 52 psf with damping; an increase 
of 150%. At Moo = 0.5, damping increased the flutter 
dynamic pressure by 50 psf from 240 psf to 290 psf, an 
increase of 21%. At Moo = 0.9, the increase in flutter 
dynamic pressure due to damping is 50 psf, an increase of 
36%. At Moo = 0.95, the increase in flutter dynamic 
pressure due to viscous damping is 62%, from 50 psf to 
81 psf. Thus the calculated flutter boundary for the AFW 
in air is sensitive to damping, varying from a moderate 
sensitivity at Moo = 0.50 to a strong sensitivity at Moo = 
0.93. 

In Fig. 12, the resultant flutter dynamic pressures at 
Moo = 0.9 and Moo = 0.93 for 0 degrees angle-of-attack and 
a viscous damping of 0.015 can also be seen. Decreasing 
the angle-of-attack from 1.5 to 0 degrees results in a 
slight stabilizing effect at Moo = 0.9 and a significant 
stabilizing effect at Moo = 0.93. An angle-of-attack 
variation performed at Moo = 0.5 revealed no difference in 
flutter dynamic pressure, as would be expected. 

Experimental Results 

During wind-tunnel testing of the AFW, a subsonic 
flutter point was encountered at Moo = 0.5 and a q=220 
psf, but it was considered to be antisymmetric and thus 
cannot be compared with the current symmetric analysis. 
It is speculated that the symmetric and antisymmetric 
flutter boundaries arc separated subsonically with the 
antisymmetric set being the most flutter critical in this 
regime; the two instabilities appear to be close, however, 
at transonic Mach numbers. 

During transonic flutter testing, three flutter points in 
the Mach number range from 0.9 to 0.93 were 


encountered. Figure 13 presents the CAP-TSD predicted 
flutter boundary at 1.5 degrees and .015 viscous damping 
(Fig. 12), the predicted linear (doublet lattice) flutter 
boundaries for symmetric and antisymmetric modes 
and the four experimental flutter points. The Moo = 0.9 
and Moo = 0.93 experimental points were also considered 
to be antisymmetric and are, again, not comparable with 
the present results. The Moo = 0.92 result, however, was 
identified as a symmetric instability, which compares 
extremely well with the CAP-TSD prediction. The 
experimental flutter frequency was about 8 Hz and the 
analytical flutter frequency was 7.8 Hz. As Fig. 13 
shows, the linear analyses predicted different trends with 
increasing Mach number. 

The no-fluttcr track in the tunnel (shown in Fig. 13) 
indicates that the bottom of the experimental transonic 
Huttcr "dip" was at about Moo = 0.93 and a dynamic 
pressure of 146 psf, much higher than that predicted by 
CAP-TSD analysis. This is not all that surprising since 
the discrepancies between TSD theory and experiment that 
exist in the heavy gas at Moo = 0,95 could be occuring at 
Moo = 0.93 in air. That is, at Moo = 0.93 in air, the flow 
could be predominantly viscous, rendering TSD theory 
inadequate at this condition. Unfortunately, there is no 
pressure data in air to verify this. This would explain 
why CAP-TSD is accurate in the region between Moo = 
0.9 and 0.92 (shocks, but possibly limited amounts of 
viscous flow) and why it is inaccurate at just a slightly 
higher Mach number, Moo = 0.93, where viscous flow 
may be dominant. A conclusive answer to this 
discrepancy at Moo = 0.93 requires additional 
investigation. 



Mach Number 


Fig. 13 Comparison of linear, nonlinear, and experi- 
mental duller boundaries in air. 
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The goal of this study was to define the Lransonic 
flutter boundary of the AFW wind-tunnel model, for 
guidance during flutter testing, and to evaluate the flutter 
prediction capability of CAP-TSD for a complete and 
realistic configuration. The static aeroclastic and dynamic 
behavior of the AFW was investigated and compared with 
experiment. 

The static acroelastic procedure developed was shown 
to be reasonably accurate. The accuracy of the procedure 
is reduced as control surfaces are deflected, probably due to 
viscous and vortex flows not addressed by TSD theory. 
As expected, the accuracy of the static aeroelastic solution 
at a given dynamic pressure and Mach number depends on 
how well the TSD assumptions represent the flow at that 
condition. 

Flutter results in the heavy gas revealed a moderate 
transonic flutter dip at Moo = 0.95, while the flutter 
boundary in air resulted in a steep dip near Moo = 0.93 due 
to the increase in mass ratio. The flutter boundary in air 
demonstrated increased sensitivity to damping and angle- 
of-attack variations at transonic Mach numbers. The 
resultant flutter boundary provided valuable guidance 
during flutter testing of the AFW, as demonstrated by the 
excellent match between theory and experiment in air at 
Moo = 0.92” The CAP-TSD results at Mpo = 0.93, 
however, were, inaccurate. At this Mach number, viscous 
effects may be dominating the flow, rendering TSD theory 
inapplicable. Until viscous acroelastic codes arc fully 
developed, however, the use of the CAP-TSD code for 
predicting the transonic aeroelastic behavior of a flexible 
wind-tunnel model, the AFW, was very valuable. 
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